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ABSTRACT 

We  determine  the  self-similar  solutions  of  the  Cauchy  problem 

vt  *  *<v>xx'  x  e  R'  1  >  0  , 
v(x,0)  -  g{x) 

for  the  constitutive  function  #(v)  *  max(0,v)  and  the  model  datum 

f  (p+)xY,  x  >  0 


-(p_)|x| x  <  0 
t  the  unique  solution  of 

"<P_HXIY»  x  <  "K^t 

x  > 

/t 


where  y»P±  >  0.  It  is  shown  that  the  unique  solution  of  (P),  (D)  is 


v(x,t) 


where 
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i/2  exp(iri  — ) 
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and  ic  is  implicitly  determined  by  the  equation  !  Avail  >13  it.y  Co 
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Dv < ♦ )  denotes  the  parabolic  cylinder  function  of  index  v.  ICI — - 1 _ 
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cylinder  functions 
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SIGNIFICANCE  AND  EXPLANATION 
Degenerate  parabolic  equations  arise  in  the  description  of  melting 
processes*  gas  dynamics  and  certain  biological  models.  The  interfaces 
corresponding  to  degeneracies  in  the  constitutive  function  usually  separate 
different  media  in  the  physical  problem. 

The  particular  problem  stated  in  the  abstract  is  related  to  nonlinear 

diffusion  equations  with  nonmonotone  constitutive  functions* as  has  been 

s - - - - - - S 

discussed  in  [HN1-3] .  in  this  report  -we  obtain  self-similar  solutions  for  (P) 
for  a  class  of  model  initial  data.  The  qualitative  behavior  of  these 


solutions*  in  particular  of  their  interfaces*  is  typical  of  the  situation  in 


more  general  problems.  In  a  subsequent  report  with  Vasques  [HNV]  we  use  such 


self-similar  solutions  as  comparison  functions  to  study  the  regularity  and  the 


behavior  for  small  time  of  the  interfaces  for  problem  (P)  with 

pi"  pH  pkt' 

0(v)  ■  raax(0,Xv) )  where  Jf  is  strictly  monotone  increasing. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


SELF  SIMITAR  SOLUTIONS  FOR  A  DEGENERATE  CAUCHY  PROBLEM 


Klaus  HSllig  and  John  A.  Nohel 


(P) 


1.  INTRODUCTION  AND  RESULT.  In  this  note  we  study  certain  aspects  of 
the  degenerate  Cauchy  problem 

vt  "  *(v)xx'  x  e  *•  1  >  0  , 
v(x,0)  -  g(x) 

for  the  constitutive  function  v)  -  max(v,0).  We  assume  that  the  initial 
data  g  are  smooth  on  R\{0}  with  at  most  polynomial  growth  at  infinity  and 
satisfy 

xg(x)  >0,  x  ft  0  , 


g(0)  -  0  . 

Problems  of  this  type  arise  as  convexifications  of  diffusion  equations 
with  nonmonotone  constitutive  functions  as  has  been  discussed  in  [HN1 ] .  The 
behavior  of  solutions  for  (P)  is  similar  to  the  one  phase  Stefan  problem 
where  g(x)  =  -1  for  x  <  0.  Existence  and  uniqueness  of  weak  solutions  of 
(P)  follows  from  nonlinear  semigroup  theory  [BCP,E] .  Moreover,  using  standard 
approximation  arguments  one  can  show  the  existence  of  a  continuous  monotone 
decreasing  free  boundary  t  ♦  s(t)  where  v(s(t)+,t)  “  0. 

The  pair  (v,s)  satisfies  the  free  boundary  problem 

vt  -  v^,  x  >  s(t),  t  >  0  , 
v(x,0)  -  g(x),  x  >  0  , 

v(s(t)+,t)  ■  0  ,  (P) 

g(s(t) )s' (t)  -  vx(s(t)+,t)  , 
s(0)  ■  0  . 


Sponsored  by  the  United  States  Army  under  Contract  No,  DAAG29-80-C-0041 . 


*  %•  *, 
-  V  * 


V. V  J*. V  -• 


Conversely,  the  solution  v  of  (P)  extended  by  v(x,t)  ■  g(x)  for  x  <  s(t) 
is  a  weak  solution  of  (P). 

We  are  interested  in  the  regularity  and  the  qualitative  behavior  of  the 
free  boundary  for  small  t.  Here  we  consider  only  the  model  datum 


P+x  * 


x  >  0  , 


^  I  Y 

[  -pJx|T,  X  <  0  , 

where  p±,Y  >  0  are  given  constants.  For  this  case  problem  (P)  has  self¬ 
similar  solutions  which  can  be  determined  explicitly.  The  problem  (P)  with 
more  general  initial  data  will  be  included  in  a  joint  paper  with  j.  Vazquez 
[HNV]  where  we  use  such  self-similar  solutions  as  comparison  functions. 

Proposition.  For  the  model  datum  (D)  problem  (P)  has  the  unique  self¬ 


similar  solution 


v(x,t)  -  tY//2i|>(— ),  x  >  s(t),  t  >  0 

✓t 


where  t|/(«)  is  the  unique  solution  of  the  ordinary  differential  equation  (1) 
below,  satisfying  the  initial  conditions  (2)  and  condition  (3)  at  infinity. 
The  free  boundary  is  s(t)  -  -x/t,  t  >  0.  For  p  :*  p+/p_  >  0  given,  x  >  0 
is  the  unique  solution  of  the  equation 


where  D_^_ ^ ( * )  is  the  parabolic  cylinder  function  of  index  (-y  -  1). 

In  the  forthcoming  paper  with  Vazquez  [HNV]  we  will  use  the  above 
Proposition  to  analyse  (P)  for  more  general  constitutive  functions  and  more 
general  data;  we  assume  there  that  $  is  smooth  on  [0,*)  with  0  < 
c  <  <  C,  $(v)  =  0  for  v  <  0  and  that  the  more  general  data  satisfy 

G(x)  -  g(x)  +  o(|x|Y)  ( | x {  ♦  0)  , 
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where  g  is  the  model  datum.  Then  it  will  be  shown  that 


s( t)  -  -<4'(0+)t  +  o(/t)  (t  \  0) 

with  te  defined  as  before.  For  $(v)  -  max(0,v)  and  Y  ■  1  we  proved  in 
[HN2,3]  the  stronger  result 

s(t)  -  -te/t  +  0(t1/2+a)  (t  V  0) 
for  any  a  <  1/2.  In  this  case  (8)  reduces  to 

2  3  2  • 

p  -  exp(-j-)  /  exp(-y2/4)dy  . 

2  4  *  -X 


Note,  that  a  first  order  expansion  of  the  equation  g(s(t))s'(t)  *  vx( s( t) ,t) 
in  (f)  formally  yields  a  different  result,  namely 

(p_)s(t)s'(t)  -  <p+)  +  •••  » 


which  yields 


s(t)  *  -  /2pt  +  •••  , 

whereas,  e.g.  for  p  -  1,  k  -  .9034  ...  *  /I.  The  reason  for  this  apparent 
inconsistency  is  that  all  derivatives  of  v  become  singular  at 
(x,t)  ■  (0,0) »  in  particular  vx  is  not  continuous  at  this  point. 


2.  PROOF  or  THE  PROPOSITION.  Substituting  v(x,t)  -  tY/2<i(x/Vt)  in  (P) 
one  sees  that  4>  aust  satisfy  the  linear  differential  equation 

2*"(£)  +  W(5)  -  Yl»(5)  -  0  for  C  >  -k  ,  (1) 

subject  to  the  initial  conditions 


Y+1 

♦  (-<)  -  0,  ♦'(-(£)  -  (pj  ^ —  (2) 

and  k  is  related  to  p^  via 

lim  rY*(5)  -  P+  .  (3) 

5-4- 


The  free  boundary  is  given  by 

s(t)  *  -ic/t  . 

Equation  (1)  can  be  solved  explicitly.  Put  x  “  £//2  and  w(£)  ■  f(x). 
Then  ( 1 )  becomes 

w"(x)  +  xw' (x)  -  Y*(x)  *  0  .  (4) 

Setting  w(x)  y(x)exp(-x2/4)  we  obtain 

2 

y"(x)  -  (j  +  y  +  f-)y(x)  -  0  .  (5) 

This  differential  equation  has  the  general  solution  [B,  p.  116-117] 
y(x)  -  ^(x)  +  b^D^ (ix)  (-•  <  x  <  •,  y  >  0)  , 

where  Dv(»)  is  the  parabolic  cylinder  function  of  index  v.  Thus  the 
general  solution  of  ( 1 )  is 


+  (C)  “ 


(6) 


for  -•  <  £  <  •  and  Y  >  0*  To  impose  the  initial  conditions  (2)  we  need  the 


formulae  (above  ref.  p.  119) 
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Then  the  initial  conditions  (2)  yield  the  pair  of  equations 


V-y-.t-  %  *  "A1’  %  * 0 

•biDr(-  7iJ  +  iTbA-i('  t=>  -  £Jj-  •*»(?->  • 


(7) 


Y+l 

K  f  1C 


Because  (5)  is  of  self-adjoint  form  the  Wronskian  of  D^(»),  )  is 

constant , 

1+1 


W(D 


-y-1<,)»Dy(*)>  S  •*«Xp[rT“)»i]  * 


Thus 


b1  (IC) 


b2(*c) 


(P_)kY+1D  (-  ■“) exp(~— ) 
_ 1  12 

i/2  exp[(-I-~-4ri] 


(p_)icY+1D_Y^1l-  ”) sxpl'g— ) 
i/T  exp[(^~-^)iri] 


(8) 


and  (6)  with  b^,b2  given  by  (8)  is  the  solution  of  (1)  satisfying  the 
initial  conditions  (2).  To  compute  the  limit  in  (3)  we  use  (see  above  ref.  p. 
122) 
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Dv(«)  -  c  exp(-  |-)[1  +  0(  |  z  | "”2 ) ]  aa  |z|  ♦  •  , 
which  is  valid  for  -  —•  <  arg  z  <  .  Thus  for  5  e  R,  y  >  0, 

2  -y- 1 

D-y-i C-—3  -  «p(-  |-)(~)  n  +  o(|er2)], 

|DYl^|  -  •xp(|r)t”z)  |«p(t1)|m  +  o(|e!‘2)],  5  - 

Substitution  of  (10)  and  (8)  into  ths  general  solution  (6)  yields 


y 

♦(5)  «  bjOOexp^!)^:)  [1  +  0(D),  C  ♦  . 

r2 


From  formula  (8)  we  see  that 


b200exp(±IH-)  -  ~  »cY+1D_y_1(-  ~)exp(~) 


Imposing  the  asymptotic  condition  (3)  and  using  (11),  (12)  we  finally  obtain 


lim  H- 

5 


which  yields  the  equation  (B). 

To  complete  the  proof  of  the  Proposition  we  have  to  show  that  given  any 
p  >  0,  (B)  is  uniquely  solvable  for  k.  From  [BO,  p.  573] 


p(<)  a  TTCT- yT  ^  exp(^H 


22r(^) 


n-0  (2n  +  1)1  VT 


is  an  analytic  function  of  it,  -*•<*<•#  Y>-1>ao“al”*' 

*n+2  *  (’f  +  1)*„  *  f  (»  '  “1  D-T-,<0)  *  2(y+1)/2’(i  +‘  ij  '  S1”°' 

the  coefficients  a  are  positive,  D  is  a  positive,  strictly 

T  /2 

increasing  function  of  tc  for  0  <  ic  <  ” ,  and  by  (14)  so  is  p(ic).  Moreover 
p(0 )  *  0.  Also  using  (14)  in  (E) 


2  2  r(i  +1) 


( K  ♦  0+) 


Moreover,  [BO,  p.  574] 


Y  2 

^2*  ,K  \T _ { K  y 


and  therefore,  from  (E) 
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x  >  0 


x  <  0 


where  Y ,p  >  0.  It  is  shown  that  the  unique  solution  of  (P) ,  (D)  is 


v(x,t) 


“(P_) |x|Y,  x  <  -e/t 


tT/2»(^=J,  » >  -c/t 
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where 


2 

♦  (5)  -  1»>1(x)D.y.1(Jz)  +  b2{l<)DY(“|)]*xPl-  |“J  * 

P.<Y+1Dy(-  ^|jexp(|-j 


b1  (k  ) 


i/T  exp(^-  -~-1  ti) 


p_<’ 


«xp(|-j 


b  (< )  "  -  v  .  , 

i/T  exp(*i  — ) 


and  ic  is  implicitly  determined  by  the  equation 

Y+1 


P  P+/P-  -  (— ]  d_y_-|1"  ~r)exptg-j  » 


D  (•)  denotes  the  parabolic  cylinder  function  of  index  v. 


